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ABSTRACT 

The Planck satellite has recently completed an all-sky galaxy cluster survey exploiting the thermal 
Sunyaev-Zekdovich (SZ) effect to locate some of the most massive systems observable. With a median 
redshift of (z) = 0.22, the clusters found by Planck at z > 0.3 are proving to be exceptionally massive 
and/or disturbed systems. One notable Planck discovery at z = 0.645, PLCK G147.3-16.6, has an 
elongated core and hosts a radio halo, indicating it is likely in the process of merging. We present 
a 16".5 resolution SZ observation of this high .-z merger using the Goddard-IRAM Superconducting 
2-Millimeter Observer (GISMO), and compare it to X-ray follow-up observations with XMM -Newton. 

We find the SZ pressure substructure is offset from the core components seen in X-ray. We interpret 
this as possible line of sight temperature or density substructure due to the on-going merger. 

Keywords: Cosmology: observations — galaxies: clusters: individual (PLCK G147.3-16.6) — galax¬ 
ies: clusters: intracluster medium — Galaxies: clusters: general — cosmic background 
radiation — X-rays: galaxies: clusters 


1. INTRODUCTION 

Forming from the largest fluctuations in the primor¬ 
dial matter power spectrum, galaxy clusters are among 
the most massive gravitationally-bound objects. There¬ 
fore, the distribution of clusters as a function of mass and 
redshift provides sensitive cosmological probes. Surveys 
spanning the electromagnetic spectrum are planned or 
underway to catalog clusters across their formation his¬ 
tory. Recent efforts exploiting the redshift-independent 
surface brightness of the Sunyaev-Zekdovich effect (SZ; 
Sunyaev & Zekdovich 1972) in particular have detected 
~1000 previously-unknown clusters (see Fowler et al. 
2010; Carlstrom et al. 2011; Planck Collaboration I 2011; 
Planck Collaboration VIII 2011; Hasselfield et al. 2013; 
Planck Collaboration XXIX 2014; Bleem et al. 2015). 

The Planck satellite has completed the first all-sky 
cluster survey since ROSAT (see e.g. Romer et al. 1994; 
Voges et al. 1999; Bohringer et al. 2000). Planck , how¬ 
ever, is not well-suited for the discovery of high-z sys¬ 
tems, whose arcminute-scale SZ signals are heavily di¬ 
luted inside Planck’s 7/3—9/7 beams at the detecting 
2 & 3 mm bands. As such, Planck detects only the 
most prominent, rare systems at high- z. The Planck 
XMM-Newton cluster validation program (Planck Col¬ 
laboration IV 2013) used the 15.5-month nominal survey 
data to identify likely cluster candidates and understand 
Planck s selection function. It suggests that the high- z 
detections are likely dynamically-disturbed massive sys¬ 
tems, which are far from being virialized and, on average, 
less X-ray luminous than X-ray selected clusters of the 
same mass. 
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Here we report high-significance 16".b resolution SZ ob¬ 
servations of a disturbed cluster from the final Planck 
XMM-Newton cluster validation program, imaging it 
with nearly 20 x better resolution than its original un¬ 
resolved detection. These new data, from the Goddard- 
IRAM 2-Millimeter Observer (GISMO; Staguhn et al. 
2006) on the 30-meter Institut de Radioastronomie Mil- 
limetrique (IRAM) Telescope 1 on Pico Veleta, Spain, re¬ 
veal complex pressure substructure in this merging sys¬ 
tem and underscore the power of subarcminute SZ follow¬ 
up. 

We summarize the known cluster properties in Sec¬ 
tion 2, discuss the new observations in Section 3, and 
present the results of our analysis in Section 4. We 
adopt a ACDM cosmology with 9^ = 0.3, = 0.7, 

and H 0 = 70 km s -1 Mpc -1 throughout this paper. At 
the redshift of PLCK G147.3-16.6 (z = 0.645; Planck 
Collaboration IV 2013), 1" corresponds to 6.9 kpc. 

2. PLCK G147.3-16.6 

PLCK G147.3-16.6 is a massive cluster at z=0.645, 
discovered at a signal to noise ratio S/N=AA1 in the 
nominal 15.5-month Planck mission. X-ray follow-up ob¬ 
servations in the Planck XMM -Newton validation pro¬ 
gram (Planck Collaboration IV 2013) reveal an ex¬ 
tended, double core morphology, while optical observa¬ 
tions with Gemini show no cD galaxy dominating the 
cluster field. More recently, 610 MHz observations with 
the Giant Metrewave Radio Telescope (GMRT) found 
that PLCK G147.3-16.6 hosts a 0.9 Mpc radio halo (van 

1 This work is based on observations carried out with the 30- 
meter IRAM Telescope. IRAM is supported by INSU/CNRS 
(France), MPG (Germany) and IGN (Spain). 
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Table 1 

Observations 


Observatory 

Date 

Project Code 

Pointing (J2000) 

R.A. Dec. 

Clean Exposure Time 
(ksec) 

GISMO 

GMRT 

XMM-Newton 

XMM-Newton 

7-9 Apr 2014 
18 Jan 2013 

27 Aug 2012 
22 Aug 2011 

235-13 

23.013 

0693661601 

0679181301 

02:56:20.0 

02:56:25.2 

02:56:23.8 

02:56:25.3 

+40:17:21.0 

+40:17:18.7 

+40:17:28.0 

+40:17:18.7 

16.9 

22 

41.7/42.1/32.9 a 

15.4/15.7/8.6 a 


Note. — Aim points and unflagged exposure times for the observations of PLCK G147.3-16.6 
included here. 

a Times for the MOS1, MOS2, and PN detectors of XMM-Newton EPIC, respectively. 


Weeren et al. 2014), placing it among the highest redshift 
radio halos known. The (re-)acceleration processes that 
create radio halos are thought to occur predominantly in 
mergers during and after the first core passage (Feretti 
et al. 2012; Brunetti & Jones 2014), although at least 
one relaxed, cool-core cluster is known to host a radio 
halo (Bonafede et al. 2014). The disturbed X-ray mor¬ 
phology and lack of a dominant cD galaxy indicate that 
PLCK G147.3-16.6 belongs to the traditional category of 
merging clusters hosting radio halos. 

3. OBSERVATIONS 

In this section we describe the resolved GISMO SZ ob¬ 
servations and new XMM-Newton data presented in this 
work. Table 1 summarizes the observations presented 
here, including those from the GMRT reported in van 
Weeren et al. (2014) and the previous XMM-Newton ob¬ 
servations in Planck Collaboration IV (2013). 

3.1. GISMO 

We observed PLCK G147.3-16.6 for a total of 4.7 hours 
(16.9 ksec) in April 2014 using GISMO, an 8x16 el¬ 
ement array of transition edge sensors. From the 30- 
meter IRAM Telescope, GISMO provides a 1/8 x 3/8 in¬ 
stantaneous field of view with 16"5 resolution at 150 GHz 
(2 mm). The GISMO data were reduced with CRUSH 2 
(ver. 2.22-1; Kovacs 2008), which was optimized to re¬ 
cover extended, diffuse signals from the atmosphere- 
dominated bolometer data. 

We estimate a total calibration uncertainty ^7-9% due 
to the imperfect knowledge of the line of sight opacities, 
based on repeated observations of Mars and Uranus (see 
Staguhn et al. 2014, for a detailed description of the ab¬ 
solute calibration). 

We deconvolved the resulting image with the measured 
point-source response of the reduction process; i.e., we di¬ 
vide the Fourier transform of the map by the 2-D transfer 
function, and back-transform. The resulting deconvolved 
image is shown in Figure 1 (left). 

The cluster field was observed using a combination of 
3-5' alt-azimuthal Lissajous patterns, yielding a median 
noise ~0.3 mJy/bm within the central 4' diameter, and 
coverage extending to approximately 6f3x8f3 area over¬ 
all. 

The noise in each map pixel was propagated from noise 
measured in the residual detector timestreams. Non¬ 
white (covariant) features, such as residual 1//, is spa¬ 
tially invariant, and hence fully captured by an appropri¬ 
ate noise re-scaling, which we determined by the ratio of 

2 http://www.submm.caltech.edu/~share/crush 
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measured-to-expected deviation, ((Sij / a ij) 2 ) > 0ut ' 

side of the approximate cluster center (at r>l/5). The 
resulting noise map, shown in Figure 1 (middle), provides 
a fair measure of the true map noise for GISMO, with 
no apparent transient noise. We find the peak SZ decre¬ 
ment at >10-<7 significance, and the overall detection is 
significant at >3-<j in every beam within the central > 2! 
of the map. 

The two-dimensional transfer function of our GISMO 
data, shown in Figure 2, was obtained by inserting a faint 
point-like test source into jackknife realizations, which 
are reduced the same way as the cluster. Thus, we ensure 
that the test dataset has the same noise properties as the 
actual dataset, and therefore the test source undergoes 
the same filtering steps as our cluster, even with adap¬ 
tive pipeline steps such as noise whitening. We averaged 
the response over 100 jackknife realizations to suppress 
the low-level sky-noise present in the individual realiza¬ 
tions. The transfer function, obtained as the ratio of the 
observed 2D spatial spectrum of the response to the un¬ 
derlying spectrum of the test source, characterizes the 
pipeline’s response to arbirtrary structures. 

Our transfer function is not circularly symmetric due 
to common mode subtraction along various correlated 
detector groups. Because CRUSH normalizes maps to 
preserve point-source peak fluxes (i.e. to keep the point 
source weighted mean response unity by definition), the 
response in the raw GISMO map diminishes gradu¬ 
ally from above unity at the short spatial scales (20"; 
3 aremin -1 ), to <1 at scales >2', which is comparable 
to the instantaneous field of view. The azimuthally- 
averaged transfer function is also shown in Figure 2 
(lower panel). 

The raw map deconvolved by the transfer function pro¬ 
vides an accurate representation of the underlying 2-mm 
flux distribution of PLCK G147.3-16.6 up to the 5' scales 
shown. The reduction and deconvolution algorithm was 
tested on both a simulated point source and a simulated 
cluster model, and it accurately reproduced the expected 
fluxes and profiles for both (see Fig. 4). The zero level 
of the deconvolved map is estimated using the mean flux 
level outside the cluster decrement (at radii R > 1.5'; see 
Section 4). 

3.2. XMM-Newton 

The discovery of PLCK G147.3-16.6 was confirmed 
through a 16.9 ksec X-ray observation with the XMM- 
Newton European Photon Imaging Camera (EPIC) in 
August 2011, as part of the Planck XMM-Newton valida¬ 
tion program (Planck Collaboration IV 2013). An addi- 
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Figure 1. Left: Deconvolved GISMO flux map, smoothed by a 12 // Gaussian to a resolution of 19. ,, 2 (depicted in the lower left corner), 
showing the SZ decrement (mjy/bm) toward PLCK G147.3-16.6. Middle: GISMO noise map with nearly uniform coverage in the 
region where the cluster is best detected. Right: Signal-to-noise (S/N) map (the ratio of left and middle images), with contours at 
S/N = [4, 2,—2,—4,—6,—8,—10] overlaid in black. X-ray contours from Figure 3 (left panel) are overlaid in white, starting at 3-cr and 
spaced at 3-cr intervals. 


GISMO Transfer Function 



Mode (arcmin 5 


GISMO Radial Transfer Function 



Figure 2. Upper: 2D transfer function used for deconvolution 
of the PLCK G147.3-16.6 GISMO observation. Note: the corner 
values at radii > 3.2 arcmin -1 are indeterminate. Lower: Radial 
profile of the above transfer function. 

tional 43.9 ksec XMM-Newton observation was obtained 
in August 2012. We calibrated both datasets using the 
Science Analysis System (SAS, ver. 13.5.0) and the most 
recent calibration files as of July 2014. The calibrated, 
cleaned event files discard periods of high intensity due 
to background particle flares. Effective exposure times 


and other observation details are summarized in Table 1. 

The images and spectra were cleaned of point sources. 
CCD4 of the MOS1 detector was operating in an anoma¬ 
lous state during observation 0679181301, and therefore 
excluded from further analysis. The merged, exposure- 
corrected X-ray surface brightness (S K ) image in Figure 3 
(left) shows that the cluster has a disturbed morphol¬ 
ogy, with an elongated core angled ~30° counterclock¬ 
wise from E-W. 

We extracted spectra for the temperature analysis us¬ 
ing the SAS tools mos-spectra and pn-spectra. Extracted 
MOS1+MOS2 spectra were co-added, and the MOS 
and PN spectra were jointly fit using xspec using cstat 
statistics. Each region contained >2000 background- 
subtracted (i.e. source) counts. With the column den¬ 
sity of hydrogen held fixed at the Galactic value, we fit 
the APEC plasma model to find the temperature of each 
region, marginalizing over abundance. Further details of 
the data reduction and analysis such as the treatments of 
the local, cosmic, and particle backgrounds are discussed 
in Bulbul et al. (2012). 

We used the XMM-Newton spectroscopic data to pro¬ 
duce a temperature ( k B T e ) map of the cluster using the 
contour binning algorithm of Sanders (2006), contbin , 
which selects regions of similar *S X above a user-specified 
S/N threshold. We used S/N > 30 yielding >2000 
source counts per region. The resulting temperature map 
is shown in Figure 3 (middle panel), with GISMO con¬ 
tours overlaid for comparison. A pseudo-pressure map 
( k B T e x a/Sx) 1S shown in Figure 3 (right). 

4. RESULTS 

In this section we compare the surface brightness of the 
GISMO SZ map with the model fit to the Planck data 
and to the properties inferred from the XMM -Newton 
observations. The thermal SZ effect traces the line-of- 
sight integral of thermal electron pressure P e . Its surface 
brightness is proportional to the Compton y parameter, 

y = [n e k B T e dl = [P e d£ (1) 

m e c z J m e c J 

where cr T is the Thomson cross-section, k B is the Boltz¬ 
mann constant, m e c 2 is the electron rest energy, n e is the 
electron number density, k B T e is the electron tempera- 
























4 


Mroczkowski, Kovacs, Bulbul et al. 



Figure 3. Left: Merged, background-subtracted XMM -Newton Sx image overlaid with GMRT 610 MHz contours (cyan) to match those 
in van Weeren et al. (2014) and with GISMO contours at S/N = [—4,—6,—8,—10] from Figure 1 overlaid (white). The central GMRT 
contours show the location of the radio halo, while those to the north of the cluster are due two to unresolved compact radio sources. The 
X-ray image, shown on a square root scale (oc density), is binned into 2. // 5x2. // 5 pixels and smoothed with a ICG FWHM Gaussian. Middle: 
XMM-Newton temperature map (in keV) using the contour binning algorithm of Sanders (2006), described in Section 3.2. Temperature 
labels show 1-cr confidence intervals in parentheses. The two white regions in the map are masked due to X-ray point sources. GISMO 
contours are depicted as in the left panel. Right: Pseudo-pressure map derived by multiplying the image (left) by the temperature 

map (middle panel) and smoothing to the resolution of GISMO. 


ture, i is the line of sight path through the cluster, and 
P e = n e k B T e . 

Planck Collaboration IV (2013) report a spherically- 
integrated Compton Y 5 oo, S ph = (5.2d=1.7) x 10 -4 arcmin 2 , 
where 

G pR500 

V 50 o, S ph = -A —2 / P e (r)4nr 2 dr. (2) 

m e c D a Jo 

Here D A is the angular diameter distance to the cluster, 
and R 500 , which for this cluster = 1042 kpc (2(5 on the 
sky), is the radius within which the average density is 
500 x greater than the critical density of the Universe at 
that redshift. Planck Collaboration IV (2013) assume the 
spherically symmetric ‘universal pressure profile’ (UPP) 
of Arnaud et al. (2010) for their model fit, shown in blue 
on Figure 4. For simplicity, we plotted only the median 
UPP fit to the Planck data, noting the error bars on 
the SZ surface brightness profile are « 13% at R$ oo? and 
~ 7% at the peak. 

The source of the discrepancy between the measured 
and modelled profiles is unknown, but may be an in¬ 
dication that the UPP is a poor fit to this disturbed 
cluster. We note that the UPP as fit by the Planck is 
treated as a matched filter function of a single param¬ 
eter, namely the mass within R$oo as computed using 
the scalings reported in Arnaud et al. (2010). Further¬ 
more, the Planck measurement of this cluster is entirely 
unresolved within R 5 oo, so the profile shown is an inter¬ 
polation of the UPP that results in the integrated signal 
measurement by Planck on scales > > R 500 . 

The GISMO S/N map shows broad qualitative agree¬ 
ment with the X-ray imaging (Figure 3, left). In the core 
however, we find the SZ and X-ray peaks are offset. An 
attempt to realign the peaks would bring the agreement 
on arcminute scales into tension. 

The observed peak offsets between the GISMO and 
XMM-Newton maps is not expected to be due to point¬ 
ing errors. XMM -Newton imaging is precise to the sub- 
arcsecond level, and the positions of bright X-ray point 
sources agree with the locations of bright WISE counter¬ 
parts. The GISMO pointing model is checked against 466 



Figure 4. SZ surface brightness profiles from the raw (dotted) 
and deconvolved (solid black) GISMO maps, with measurement 
uncertainties (dark grey) and systematic zero-level estimation un¬ 
certainty (light grey) ranges, compared with that computed from 
the median UPP (blue) fit to the Planck data (Planck Collabora¬ 
tion IV 2013). We also show the profile we recover (red with 2-cr 
uncertainties) when we insert the UPP into the jackknifed GISMO 
timestreams and analyze it the same way as the actual cluster ob¬ 
servation. 

measurements of bright point sources, observed hourly 
through several days. All are within an RMS deviation 
< 4"5. The GISMO observations combine 32 observa¬ 
tion blocks, each bracketed by independent pointings, 
thus the statistical 1-cr pointing error of the composite 
map is estimated at 0'.'8. We also note that the positions 
of dozens of bright compact sources observed by GISMO 
observations during the same observation period are re¬ 
produced within the expected accuracy. We therefore 
consider the discrepancy to be of astrophysical nature. 

We do not expect to detect any contamination by the 
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Cosmic Infrared Background (CIB). The unresolved part 
of the CIB is removed by flux zeroing outside of the clus¬ 
ter, so our fluxes are effectively referenced against the 
mean CIB level. Resolved CIB sources may be present 
in the map, but are unlikely to be detectable. A deep- 
field study with GISMO by Staguhn et al. (2014) finds no 
sources brighter than ~1 mJy (>3 -a in our map) in a sim¬ 
ilar area to our field, and put the 2 mm confusion noise at 
<50/iJy, i.e. several times below the RMS in the observa¬ 
tions presented here. The combined CMB+CIB at scales 
from 1-4' has also been measured at 150 GHz by the SPT 
to be <100/iK (George et al. 2015), corresponding to a 
signal <407/iJy/bm in the GISMO data, comparable to 
the noise level in our map. We also note that recent stud¬ 
ies by Sayers et al. (2013) and Gralla et al. (2014) both 
found the radio and submillimeter point source contribu¬ 
tion to be minimal near 150 GHz. Therefore, CIB/CMB 
contamination in our map is expected to be small. 

SZ and X-ray imaging are sensitive to the line of sight 
integrals of pressure (Eq. 1) and density squared, re¬ 
spectively; for bremsstrahlung emission, X-ray surface 
brightness has only a weak dependence on temperature, 
5 X oc fn e 2 Differences between the SZ and X- 

ray maps (y and 5 X ) can therefore be due to temperature 
substructure or to the differing line of sight distribution 
of the gas (e.g. clumping or asphericity). 

The location of the SZ signal at >6-cr broadly agrees 
with the location of the hottest gas found in the temper¬ 
ature and pseudo-pressure maps (Figure 3, middle and 
right, respectively), but we find no clear evidence for 
shock-heated gas at the resolution of the X-ray spec¬ 
troscopy. For the high-significance (>6-cr) SZ region, 
we extracted spectra from both XMM-Newton observa¬ 
tions, and X-ray counts within R 50 o, using the same fit¬ 
ting procedure and plasma model as for deriving the 
temperature. We find only a marginal enhancement 
of T x = 11.331^;^ keV over the global temperature 
^x, 5 oo = 8.74 ^o'. 56 keV, which agrees with that found 
by Planck Collaboration IV (2013). 

This leaves the possibility that the SZ/X-ray offset is 
due to an irregular gas distribution along the line of sight 
or the breakdown of the assumption that pseudo-pressure 
and SZ features should directly match. The so-called 
‘slab approximation’, which treats the line of sight tem¬ 
perature as isothermal in each spectroscopic bin and as¬ 
sumes the path length through the cluster is a constant 
(e.g Mroczkowski et al. 2012; Planck Collaboration X 
2013), may not hold for complicated merger geometries. 

5. CONCLUSIONS 

We present the first high significance maps of the SZ 
effect with GISMO, revealing substructure in the Planck- 
selected cluster PLCK G147.3-16.6. The core morphol¬ 
ogy mimics the appearance of the X-ray observation (re¬ 
ported here and in Planck Collaboration IV 2013), but is 
notably offset from their X-ray counterparts. The pres¬ 
ence of the giant radio halo reported in van Weeren et al. 
2014 further supports the hypothesis that this system is 
likely a merger. 

This GISMO observation demonstrates that a compa¬ 
rable level of detail in a moderately high-z cluster can 
now be obtained from large, ground-based telescopes in 
a similar amount of time as that currently required for 


X-ray observations. This adds to a small but growing 
number of instruments that have imaged the SZ effect 
at resolutions better than 20", which include such in¬ 
struments as Nobeyama (e.g. Komatsu et al. 2001; Ki- 
tayama et al. 2004), MUSTANG (e.g. Mason et al. 2010), 
CARMA (e.g. Plagge et al. 2010), and NIKA (e.g. Adam 
et al. 2014, 2015). New and future observations with 
GISMO and GISMO-2 (Staguhn et al. 2012) will probe 
high -z cluster mergers further, or confirm cluster candi¬ 
dates from SZ, X-ray, and optical surveys. 
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